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Abstract 

An extension of Hofava-Lifshitz gravity was recently proposed in order to address 
the pathological behaviour of the scalar mode all previous versions of the theory 
exhibit. We show that even in this new extension the strong coupling persists, 
casting doubts on whether such a model can constitute an interesting alternative to 
general relativity (GR). 
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A new perspective towards addressing the quantum gravity puzzle has received 
increased attention recently: Hofava proposed a theory of gravity which does 
not respect Lorentz invariance [1]. Instead, it assumes the existence of a pre- 
ferred foliation by 3-dimensional constant time hypersurfaces, which splits 
spacetime into space and time. This allows to add higher order spatial deriva- 
tives of the metric to the action, without introducing higher order time deriva- 
tives. This is supposed to improve the ultraviolet (UV) behavior of the gravi- 
ton propagator and render the theory power-counting renormalizable without 
introducing ghost modes, which are common when adding higher order cur- 
vature invariants to the action in a covariant manner [2]. The new theory is 
supposed to be an adequate UV completion of GR, above its natural cutoff 
scale, which is of the order of the Planck scale. 

Such a theory cannot be invariant under the full set of diffeomorphisms, but 
it can still be invariant under the more limiting foliation-preserving diffeomor- 
phisms, t — > t(t), x l — > x l (t, x l ). In this setting, it is natural to consider the 
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Arnowitt-Deser-Misner (ADM) decomposition of spacetime 

ds 2 = -N 2 c 2 dt 2 + gij {dx l + rfdtXdx? + N'dt). (1) 
Defining the extrinsic curvature as 

K« ^ {/A, -- V,-Vr- V,,V, ) . (2) 

the action of the theory is of the form 

M 2 r 

S=^ J d 3 xdtN^/g[K ij K ij - XK 2 - f) , (3) 

where M p i is the Planck mass, g is the determinant of the spatial metric gij 
and A is a dimensionless running coupling. "V depends only on and its 
spatial derivatives. 

Power-counting renormalizability requires that "V include at least sixth order 
terms in derivatives of g^. To see this, it is better to switch from the c = 1 
units used above, to units that impose the scaling 

[dt] = [«r 3 , [dx] = [K]-\ (4) 

where k is a placeholder symbol with dimensions of momentum. This makes 
the couplings of the kinetic term and the sixth order derivatives dimensionless 
and renders the theory power counting renormalizableQ Switching back to 
c = 1 units, which are more suitable for discussing the infrared behavior, "V 
has the general structure 

f = 2A- R + M-?0(R 2 ) + M-fO(R 3 ), (5) 

where, R is the Ricci scalar of g^ and 0(R 2 ) and 0(R 3 ) denote all terms 
that can be constructed from g^ and its spatial derivatives up to dimensions 
4 and 6 respectively. A plays the role of a cosmological constant and can be 
easily expressed in terms of a dimensionless coupling as A = goM p \/2, whereas 
the coefficient of R has been set to —1 by a suitable rescaling of the coordi- 
nates. Note that the natural scale suppressing higher order operators here is 
the Planck scale M p \; we will return to this point later. From a phenomeno- 
logical perspective, the hope is that A will flow to 1 in the IR, which is the 
general relativistic value, and considering that the higher order terms will be 
suppressed, Lorentz invariance will be restored as an emergent symmetry. 

In principle "V should include all possible terms up to a suitable dimension. 
Hof ava proposed a specific structure for "V inspired by condensed matter sys- 
tems, dubbed "detailed balance", in order to reduce the number of terms [1]. 

1 See also [3] for a more extensive discussion of this principle, using the simpler 
case of a scalar field as an example. 
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Since its physical motivation is unclear and it leads to the wrong sign of the 
cosmological constant (see however [4]), we will not restrict ourselves to it 
here. In any case, our results will be independent of the specific choice of "V . 

Another restricted version of the theory is the one based on the assumption of 
"projectability", N = N(t). This was proposed as a way to match the reduced 
symmetry and still be able to set N — > 1 using gauge transformations. This 
assumption significantly reduces the number of terms in "V . It has been used 
in Refs. [5,6] as a simplicity assumption in order to construct and study the 
most general action in this framework I" 2 ! 

As already pointed out in Ref. [1], the breaking of general covariance generi- 
cally introduces a scalar degree of freedom. The behavior of this scalar creates 
a serious viability issue for all versions of the theory. To be specific, in the 
version with "projectability" and for maximally symmetric backgrounds, the 
scalar is (classically) unstable for A>lorA<l/3 and it becomes a ghost 
(quantum-mechanically unstable) for 1/3 < A < 1 [6,8,9]. For the same back- 
grounds, when A — > 1 the cubic interactions blow up and the scalar becomes 
strongly coupled [10,8]. 

The non-projectable version is also burdened with similar problems. Even 
though the scalar mode appears to be absent (frozen) from the spectrum for 
backgrounds that are either time-independent or spatially homogeneous (e.g. 
maximally symmetric or cosmological backgrounds [13,10]), it reveals itself for 
more complicated backgrounds [10] and signals the breakdown of the theory 
at unacceptably long distances, leading to strong coupling and instabilities. 
See also [11,12] regarding pathologies of the non-projectable case. 

Having in mind that the problems in the non-projectable case stem from the 
anomalous structure of the quadratic action for the scalar degree of freedom 
around smooth backgrounds, Bias, Pujolas and Sibiryakov recently proposed 
an extension of the theory which leads to a healthy quadratic action [14]. Our 
purpose here is to examine whether this extended version of Hofava-Lifshitz 
gravity also avoids the strong coupling issues burdening its predecessors. 

We start by reviewing the extended theory proposed in [14]. The key idea is 
to add to "V terms with up to six spatial derivatives, constructed with the 
3- vector 

a, = diN/N, (6) 

and its contractions with 3-dimensional curvature invariants of g^. Such terms 
are manifestly invariant under the symmetry of the initial theory, t — > t(t), x 1 — > 
x l (t,x l ). In principle, one can straightforwardly write down all of the new 
operators of dimensions 4 and 6 that involve a, and its contractions with 



2 See also [7] for a more cosmologically oriented study. 
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3-curvature invariants and modify the theory in the UV. These terms are nu- 
merous and we will avoid listing them here. However, note that there are more 
termj 3 ! than those contributing to the quadratic action listed in Ref. [14]. 

In the following, we will focus on the IR limit of the theory. This means that we 
will take into account operators in of lowest order, i.e. of at most dimension 
2. Apart from the first two terms in (5), there is only one more term involving 
a 1 that has to be added, 

V aid 1 , (7) 

where rj is a coupling constant, which is in principle also running, just as A. 
In fact, the full IR action, in which we will focus from now on, is (setting 
M pl = 1) 

S = ^J d 3 xdtNy/g [K^K^ -XK 2 + R + r] a^) . (8) 

To study the interactions of the theory, we will adopt a gauge similar to the 
one used in Ref. [8] and write down the scalar perturbations of the metric as 

N = e a ^ x \ N, = diP(t, x), 9lj = e 2C ^%. (9) 

The above form differs from the most general scalar perturbation possible, by 
the perturbation of gij of the form 2didjE(t, x). Such a perturbation transforms 
under the foliated diffeomorphism group as E — > E — \ [15], where x % — > 
x l + x % ■ It is then straightforward to use the freedom in x{t, x) in order to set 
E = 0. 

Substituting eqs. (9) into eq. (8) leads to the quadratic Lagrangian, which 
after partial integrations reads 

C 2 = ^(l-3A)C 2 + (9C) 2 + ^(l-A)(A/3) 2 -(l-3A)CA/3-2aAC-|aAa. (10) 

This is the quadratic Lagrangian of Ref. [8], where the projectable version of 
Hofava-Lifshitz gravity was considered, with two important differences. Since 
iV is now spacetime dependent, a is spacetime dependent as well and, there- 
fore, the penultimate term in eq. (10) is not a total derivative. Secondly, the 
extended theory has the extra coupling (7) which contributes to the quadratic 
action. 

Varying the quadratic action with respect to (5 and a, we obtain the momen- 
tum and the Hamiltonian constraints (assuming regular boundary conditions) 

A/3 = -ic, « = --C (11) 

C\ 7] 

3 Since ./V = N(t,x l ), integration by parts does not allow to drop as many terms 
as one would expect. For instance, V 2 R is no longer a surface term. 
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In the above, we used the quantity 



1-A 

whose physical significance we will see shortly 



(12) 



Substituting these constraints to the quadratic action and with appropriate 
partial integrations, we obtain the following action for the physical mode £ 



d 3 xdt 



\( 2 - v —(d0 2 



V 



(13) 



This is in agreement with Ref. [14], if we ignore the higher order operators 
(of dimension 4 and 6) which were kept there. In the original theory, without 
the r\ operator, it is evident from (13) that would be the speed of sound 
for the perturbation £. Also, from the same action, we see that without the r\ 
contribution, £ would either be unstable when < 0, or it would be a ghost 
when c 2 > 0, as mentioned above. The improved theory, however, can cure 
this problem for 

c\ < and < rj < 2. (14) 

We now turn our attention to the next order in perturbation analysis in order 
to check whether strong coupling can be avoided. The cubic Lagrangian reads 



C 3 = ^(1 - 3A)CC 2 - C(<9C) 2 - C 2 AC - (1 - 3A)CCA/3 - (1 - 3\)(d k (d k (3 

+(1 - \)A!3d k (d k !3 - IcdidjPtf&P - ■><);<), ur + ^C(A/3) 2 

--(1 - 3A)«C 2 + (1 - 3A)«CA/3 - l -ad t d 3 Pd l d 3 p + -a(Ap) 2 - a(d() 2 
2 2 2 



-2aCAC - a 2 AC + ^((da) 2 + ^a(da) 2 , 



V 



(15) 



Substituting the constraints (11), in the same fashion as above for the quadratic 
action, leads to the cubic action 



6', 



dtd^x 



/3 1 N 
+ 2 + = 



2 • d i ■ 
c£ A 



4c 



didj 
~A~' 



(16) 
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We note again that when if 77 — > 00, we obtain the couplings found in Ref. [8] 
(once the momentum constraint is substituted) for the projectable case PI 



After restoring M p \ and performing canonical normalization of the kinetic term 
in the quadratic action, eq. (13), as ( = \c^\(/M p i, all the terms but the first 
in the cubic action, eq. (16), scale as (|c^|M p i) _1 . Therefore there is strong 
coupling for c<^ — > (A — > 1), i.e. the £- interactions become strong for energies 
above the scale |c^|M p i. We note here that all terms that blow up in that limit 
have time derivatives of (. 

In order to see if the presence of the new coupling 77 can cure the strong 
coupling problem found above, we should examine whether its contributions 
could in general be used to cancel (with fine-tuning, since 77 is in principle 
running) the troublesome interactions. This can happen for the interactions 
in (16) that depend on 77. However, in these cases, the values of 77 required 
for the cancelation to take place lie outside the region < rj < 2, contrary 
to what is needed in order to have a healthy mode. Moreover, even if such a 
cancelation does happen, no value of rj can erase all divergent terms in the 
limit — > 0. The third term in (16) in particular is 77- independent and cannot 
be canceled for any value of 77. Therefore, the terms that blow up at cubic 
order cannot be removed from the action of the theory, not even with fine 
tuning of the new coupling 77. 

The higher order terms of dimension 4 and 6 which we have ignored, can also 
not be used to cancel out the troublesome interactions. This is because they 
generate new couplings involving only spatial derivatives of ( and a, but not 
p. Only the latter is related to time derivatives of ( whose interactions exhibit 
strong coupling. An effect the higher-dimensional operators could have is the 
following: if instead of assuming they are suppressed by powers of the M p \ 
(which is the natural scale here) we introduce a new scale Z exactly for this 
purpose — which is equivalent to tuning the dimensionless couplings — then, 
if Z < \c^\M p \, the perturbative expansion will be dominated by the higher- 
dimensional operators before strong coupling kicks injjO We will discuss this 
possibility, which clearly requires tuning, separately in what follows. 

So far we have shown that there is inevitable strong coupling as A — > 1, at 
least as long as a new scale Z is not introduced. It now remains to discuss 
how physically relevant is this strong coupling region of the theory. First of 
all, just by inspecting the IR action (8), it is straightforward to see, that 
Lorentz invariance and diffeomorphism invariance are recovered when A — ^ 1 
and 77 — > 0. Indeed, in the initial version of the non-projectable Hofava-Lifshitz 
gravity where 77 = 0, it was hoped that A would flow to 1 in the infrared in 

4 The projectable theory limit as noted in Ref. [14] is obtained for 77—7-00. 

5 This possibility was brought to our attention by the authors of Ref. [14] in private 
communication and was subsequently the subject of Ref. [16]. 
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order for Lorentz invariance to emerge. Clearly, this is not an option here, 
even if 77 would flow to for some reason, as the resulting theory would be 
strongly coupled at all scales. 

The only other alternative is to have a theory which exhibits Lorentz violations 
all the way to the IR regime, as A does not flow to 1 there (and neither does 77 
flow to 0). In general, this is probably more likely anyway, as there is no reason 
to expect that the renormalization group flow would work to the benefit of 
Lorentz invariance recovery. However, in that case, how far from A = I and 
rj = are we allowed to be, from an experimental perspective? Gravitational 
experiments, and also non-gravitational ones once matter is consistently added 
to the theory, should be able to provide tight constraints. This is an avenue 
worth exploring further. 

However, for the time being it suffices to consider for purely demonstrative 
purposes a specific constraint already examined in [17,14]: the one coming from 
the discrepancy between the effective gravitational constant as measured in 
a gravitational experiment G e s and that measured in cosmology, G C osmo- As 
pointed out in [14], when considering a static metric around a point source 
one can start from the quadratic Lagrangian (10) supplemented by the source 
term C m = T 00 SN = —ma8^(x), and straightforwardly derive the Poisson 
equation (after reworking the Hamiltonian constraint in the presence of the 
source coupling) 

AC = - m * 3(X \ . (17) 
_ 2(1-77/2) 1 J 

Simple comparison with the usual form of the Poisson equation reveals that 



G 



eff 



On the other hand, whenever N is not a function of time, a 1 vanishes, and this 
is the case in cosmology, where one assumes that the universe is homogeneous 
and isotropic at large scales. Thus, the cosmological equations will be insen- 
sitive to the new coupling coming from the a 1 terms, and the first Friedmann 
equation (including phenomenological matter) at low energies will be the usual 
one [6,7] with effective gravitational constant as measured in cosmology 

Gcosmo = 8tt(3A-1)' (19) 

The discrepancy between G e fr and G CO smo is constrained by measurements of 
the primordial abundance of He 4 [18]: 

\G cosmo /G eS - 1| < 1/8. (20) 

Combining this constraint with those in eq. (14) and the fact that we want 
Gcosmo > 0, we obtain that there is a triangular region of the (77, A) parameter 
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space which is allowed. For 77 = 0, we obtain the maximum upper bound for 

A as 

0<A-1<0.1. (21) 
Thus, if the theory is to be viable, A will have to be sufficiently close to 1. 

Continuing the previous analysis to higher orders, we anticipate that the inter- 
actions in the perturbation series will scale as (\c^\M p i)~ n ~ (y^l — A|M p i)~ n 
for increasing orders n (as was the case for the previous versions of the the- 
ory as well [10]). Based on the mild bound above, and assuming that no new 
scale Z is introduced, the strong coupling kicks in at the energy scale 0.2M p i, 
which is already an order of magnitude lower than the natural cutoff of GR. 
Note that the experimental constraint used here is not the only one A has to 
satisfy in principle. In fact, this is just an extremely mild constraint and other 
constraints coming from Lorentz violations would be much more stringent, 
requiring that A be even closer to 1. For example, in [19] the 77 = version 
of the theory was studied, yielding a bound A — 1 = 0.003 ± 0.014 from the 
running of the Planck mass. Another bound mentioned in [16] from the c*2 
PPN parameter yields A — 1 < 10~ 7 . These and similar bounds will push the 
strong coupling scale orders of magnitude lower and definitely deserves fur- 
ther investigation. In any case, the very existence of such of an effective cutoff 
would imply that this extended theory cannot be considered a UV completion 
of GR, which is what it was introduced for. That would be true irrespectively 
of whether the strong coupling energy scale is low enough to fall within the 
range in which we can test gravity theories. 

Let us summarize. We examined the extended version of Hof ava-Lifshitz grav- 
ity proposed in [14] as a way to address viability issues caused by the anoma- 
lous dynamics of the scalar degree of freedom that is present in the original 
version of the theory. Expanding around flat space, it was shown that, even 
thought the quadratic action of the theory is "healthy" , there are cubic inter- 
actions that blow up when A — > 1. This leaves no space for Lorentz symmetry 
to be recovered at low energies. Additionally, even if one is willing to abandon 
Lorentz symmetry altogether, experimental constraints imply that A has to 
be close to 1. The strong coupling scale, which acts as an effective cutoff of 
the theory, appears to be lower than the cutoff of GR as an effective theory. 
Considering that the motivation of the model is to constitute a renormalizable 
theory of gravity, this casts serious doubts on whether it can be an interesting 
alternative to GR. 

As mentioned earlier, a possible way out would be to introduce an new scale Z 
and use this scale instead of the Planck scale for suppressing the higher order 
operators, so that these operator could take over the perturbative expansion 
before strong coupling kicks in. However this definitely requires tuning and it 
would cause naturalness issues, as such a scale would need to be motivated by 
some additional physics. Note also that the closer to 1 observations required 
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A to be, the smaller the scale Z should be chosen to be as well. Trying to 
make A small in order to avoid Lorentz violations related to its present in 
the kinetic part of the action will effectively lower the scale at which Lorentz 
violations related to the higher order operators appear and vice versa. Thus, 
in this scenario A will be bound both from above and below, and it remains 
to be seen how wide is the range of values it can actually take. (As already 
mentioned, the experimental constraint used here for demonstrative purposes 
was an extremely mild one.) 

Another possible way for the theory to make sense despite the strong coupling 
could be through some analogue of the Vainshtein mechanism [20] in massive 
gravity. It is conceivable that such a non-perturbative restoration of the limit 
A — > 1 will take place in the present theory as well. Alternatively, one may 
hope that the renormalization group flow of A will conspire in order for A to be 
much larger than 1 around the Planck scale and at the same time to run very 
close to 1 at low energies. If such a tuning can be achieved then the improved 
theory could constitute an adequate UV completion of GR. 
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